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Abstract 

Atomic quantum gases in optical lattices serve as a versatile testbed for important 
concepts of modern condensed-matter physics. The availability of methods to charac- 
terize strongly correlated phases is crucial for the study of these systems. Diffraction 
techniques to reveal long-range spatial structure, which may complement in situ de- 
tection methods, have been largely unexplored. Here we experimentally demonstrate 
that Bragg diffraction of neutral atoms can be used for this purpose. Using a one- 
dimensional Bose gas as a source of matter waves, we are able to infer the spatial 
ordering and on-site localization of atoms confined to an optical lattice. We also 
study the suppression of inelastic scattering between incident matter waves and the 
lattice-trapped atoms, occurring for increased lattice depth. Furthermore, we use 
atomic de Broglie waves to detect forced antiferromagnetic ordering in an atomic 
spin mixture, demonstrating the suitability of our method for the non-destructive 
detection of spin-ordered phases in strongly correlated atomic gases. 



The diffraction of electron and neutron matter waves from crystalline structures is a stan- 
dard tool in crystallography, complementing x-ray techniques [lj. The advent of quantum 
gases in optical lattices has introduced a new class of highly controllable systems that mimic 
the physics of solids at lattice constants that are three orders of magnitude larger [2], and 
it appears natural to ask about a possible role for atomic matter-wave diffraction in the 
characterization of these systems j3j [4] . 

Several approaches for detecting the spatial structure of strongly correlated phases of 
ultracold atoms in optical lattices have so far been taken. These include the analysis of 
noise correlations in momentum space [6], accessible after release of the atoms from 
the lattice, as well as dedicated in situ detection methods based on optical imaging [7-9J 
and electron microscopy [10J. In this context, diffraction experiments have the potential 
to reveal important information on the existence of long-range order, such as spin-ordered 
phases in atomic mixtures [TTHT3] . in a non-destructive manner and with substantially lower 
experimental requirements. Here, optical and atomic matter- wave diffraction are equivalent 
in the sense that scattering of photons and atoms can be sensitive to both the external 
and internal state of atomic scatterers. However, there are certain advantages to using 
matter- wave probes. The de Broglie wavelength of an atomic probe can be tuned freely by 
controlling its velocity, thus precluding limits on spatial resolution and also providing access 
to Bragg resonances without the need of varying the angle of incidence. As in the optical 
case [TH [15] , matter- wave probes of high spectral brightness are readily available by using 
atoms from a Bose-Einstein condensate. 

In this work we study the scattering of a probe species, consisting of one-dimensional 
Bose gases, from target atoms confined to an optical lattice. For weak confinement, we 
observe free-particle-like, one-dimensional (ID) collisions [HI Ej between the two species, 
corresponding to inelastic band-structure excitations of the target by the incident probe field. 
The inelastic scattering is suppressed for more deeply confined target atoms, giving way to 
elastic scattering and Bragg diffraction, from which the underlying crystalline order can be 
inferred. We use this matter- wave scattering technique to characterize the introduction of 
forced antiferromagnetic ordering [18] in the system. 

Matter-wave probing of an atomic crystal 

Our experiments begin with a virtually pure Bose-Einstein condensate of 87 Rb atoms in the 
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|F, mp) = |1, —1) hyperfine ground state, prepared in an optical dipole trap (ODT) of nearly 
isotropic harmonic trapping frequency 2tt x 50 Hz in the transverse [x and y) and vertical (z) 
directions. Along the transverse axes, far-detuned attractive optical lattices (A_l = 1064 nm) 
are smoothly ramped up in 200 ms, resulting in an array of one-dimensional Bose gases, with 
trapping frequency uj z /2ti « 70 Hz along z. The final depths of these transverse lattices, 
s± = 40 (measured in units of the transverse recoil energy (/i/A^) 2 /2m, where m is the 
atomic mass), are sufficiently deep so as to energetically restrict all dynamics to the lowest 
mode of the transverse potential, and to suppress tunneling between the ID tubes on our 
experimental timescales. 

We next create a mixture of probe = |2, — 2)) and target (\t) = |1,0) or |1,— 1), 
see below) species through microwave manipulation of the atoms' internal hyperfine state. 
By adiabatically ramping up a magnetic field gradient, we then fully separate the probe 
and target atoms along the longitudinal (z) axis. We thereafter smoothly ramp up a state- 
selective lattice along z, formed by light of wavelength X z = 785 nm, between the Di and D 2 
lines of 87 Rb. This allows ( [T9142T] . see Methods) for a cancellation of light-shifts (and zero 
lattice depth) for the probe atoms, while the target atoms experience an attractive lattice 
along z, with period d = X z /2 and variable depth s z (measured in units of the longitudinal 
recoil energy E R = (h/ \ z ) 2 /2m). This allows for the target to be driven to a ID Mott 
insulator state [22j, while the probe species remains a superfluid ID Bose gas. Finally, as 
detailed below, we examine scattering processes between probe and target at a well-defined 
relative velocity, v re ^ relating to a probe de Broglie wavelength of A^b = h/mv re i. 

ID collisions and inelastic scattering 

We begin by studying collisions between the two species, as the target atoms become local- 
ized to a state-selective lattice of increasing depth. In the limit of zero lattice confinement 
along z, the collisions occurring between the \p) and \t) atoms are essentially free-particle 
collisions, with all scattering restricted to the ID tubes. Such binary collisions have previ- 
ously been studied with a single species as a quantum analog of a "Newton's cradle" [T7j, 
and with a two-species mixture in the context of damped spin impurity transport [23]. Here, 
we realize such collisions by accelerating the probe atoms into an initially separated sample 
of target atoms, which is itself at rest. 

As illustrated in Fig. [j] (a), we set the relative velocity of the collisions by controlling the 
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magnetic field gradient along z that initially separates the two species. The probe atoms 
experience a longitudinal trapping potential shifted by a distance Az, while the target 
atoms (|jP, mp) = |1>0), insensitive to the applied gradient) remain at the trap center. The 
magnetic field gradient is quickly extinguished, and after a quarter oscillation in the trap 
the probe atoms have accelerated to a nominal velocity of v® = u z Az upon interacting with 
the stationary (v® = 0) target atoms. At this point, we access their momentum distributions 
by turning off all confining potentials and separating them by applying a pulsed magnetic 
field gradient during time of flight (TOF), followed by absorption imaging. 

Typical TOF absorption images are displayed in Fig. [j] (b), for the case of free (along 
z) target atoms and probe atoms incident at v re \ — v® — v® — 2vr (^-lattice recoil velocity 
vr — h/2md ~ 5.8 mm/s). Here, and also for weak ^-lattice depths, two distinct velocity 
components at and v re \ can be seen for both species. As the target atoms are initially at 
rest and the accelerated probe atoms are initially at v re \ before colliding with each other, the 
scattering spectra show energy- and momentum exchanging binary collisions between probe 
and target atoms, i.e. reflections in the center-of-mass frame of atomic pairs. Whereas reflec- 
tion and transmission events are indistinguishable for colliding atoms of the same spin [T7j , a 
spin-mixture gives experimental access to the reflection probability [IB] . The atoms collide 
with a high kinetic energy that far exceeds the mean-field energy of either species, such 
that the collisions have free-particle character. Using a slightly uneven mixture of the two 
species (3:2 target to probe atoms), we find that 11 ± 1% of the target atoms and 14 ± 3% 
of the probe atoms are reflected, in fair agreement with the calculated reflection probability 
of R — 10% [E] for our system parameters. 

As mentioned above, much effort has been devoted to studying collisions of lattice-free ID 
Bose gases, both for the spin-polarized [TTj and spin-mixed cases [23]. In our system, we now 
investigate what happens when the dispersion relations are qualitatively different, with one 
of the species being subject to a longitudinal optical lattice. To discuss in simple terms our 
expectations for the case that the target atoms reside in a lattice of finite depth, we describe 
the momentum and energy exchange in the periodic zone scheme associated with the lattice. 
For zero lattice depth, the dispersion relations of target and probe atoms coincide, such that 
the momentum exchange is resonant. However, as illustrated in Fig. [I] (c), for non-zero 
lattice depth s z ^ 0, each collision demotes an incident probe atom from the first band 
to the ground band, and vice versa for the target atom involved in the collision. With an 
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energy mismatch SE that increases with s Z) the momentum-exchange becomes off-resonant, 
amounting to an inelastic band-structure excitation of the target. This process can occur as 
long as the energetic uncertainty AE ~ h/r associated with the finite interaction time r is 
larger than 5E. 

The observed changes of the target's TOF spectra with increasing lattice depth are ana- 
lyzed in Fig.[j](d). In the absence of collisions, ^-dependent, symmetric peaks due to optical 
diffraction are observed at v t = ±2^ as expected (for up to s z ~ 20). Collisions at the inci- 
dent velocity v® = 2vr give rise to an asymmetry between the two peaks, depending on the 
depth of the lattice. The observed asymmetry decays with the mismatch 5E ) in qualitative 
agreement with our expectation. We estimate the time for probe atoms to traverse half of 
the target (at which point we switch off the lattice) to be r ~ lO/xm/2^ ~ 0.9ms, giving 
an associated energy uncertainty AE ~ 0.3E R . This value is in good qualitative agreement 
with the observed \ je exponential decay constant of 0.35E R . Our observation of the sup- 
pression of elastic two-body collisions due to band-structure mismatch should be of direct 
consequence for recently proposed thermometry schemes in strongly correlated systems that 
are based upon the use of a lattice-free spectator species in contact with lattice-modulated 
atoms [ 24] . 

Elastic Bragg diffraction of matter waves 

While inelastic scattering events are seen to die off with increasing lattice confinement, we 
instead can expect to observe elastic scattering of probe matter waves from the crystal of 
target atoms. The distribution of elastically scattered probe atoms is determined by the 
static structure factor S(q) of the target (where hq is the probe momentum transfer). As 
S(q) is given by the square of the Fourier transform of the target's density, probe scattering 
can thus reveal information about the underlying spatial distribution of the target atoms. 
As a caveat, we briefly mention that the study of the dynamic structure factor through 
low-energy inelastic scattering, which could provide insight into correlations of the target 
atom array [4J , is largely precluded in our system due to dominant reflection at low energies 
(described below), as well as our use of a superfluid probe gas, which itself supports low- 
energy collective excitations. Such a study may thus be more well-suited to energy-deposition 
measurements based on optical Bragg spectroscopy [25j [26] . 

The static structure factor, for a deeply modulated target, can be approximated as S(q) = 
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I Ylj fj(Q) eiqZj \ 2 ? where the fj(q) are the amplitudes of isolated scatterers at positions Zj = 
jd, with j an integer. If one assumes them to be identical (fj(q) = /(g)), as in a Mott 
phase with uniform filling, the structure factor is a product of two terms, | ^-e iqZj \ 2 and 
l/(?)| 2 — K^o \e iqz \ ^o)| 2 , where (fo(z) describes the on-site density distribution of each 
scatterer. The first term determines the positions of Bragg diffraction resonances, and 
reveals the periodicity of the array of scatterers. The second term, the atomic form factor 
|/(g)| 2 , or equivalently the single-scatterer envelope function, reveals information of the 
on-site density distribution. 

In order to perform crystallographic measurements, we begin by introducing a controlled 
method to vary the relative velocity (de Broglie wavelength) of matter waves incident on a 
deeply-confined crystal. Such a wavelength scan is necessary to identify Bragg diffraction 
resonances, which are expected when a multiple of the probe's de Broglie wavelength coin- 
cides with the target's lattice spacing (i.e. 2d = nA^, or v re \ = ttor, with n the order of 
diffraction), and more generally it allows study of the dependence of scattering on the probe's 
de Broglie wavelength. As illustrated in Fig. [2] (a), the target atoms (here in the |1, — 1) 
state) are loaded into a very deep state-selective lattice initially at rest in the lab frame. 
This lattice is then made to move at a well-defined velocity vl = S(jj(X z /4tt) by detuning the 
relative frequency of the two interfering laser beams that create the lattice by an amount 
duo. The pinned target atoms, which are initially far-separated from the probe atoms, follow 
the moving lattice minima at a variable velocity v re \ for a total time T move , including a set 
delay time before the probe and target make contact. The large initial separation helps 
to avoid the formation of a quantum emulsion [27J, resulting when homogeneous mixtures 
are loaded into state-dependent lattice potentials, which can affect interspecies transport 
properties [211 EH] • 

In Fig. [2] (b), we show TOF absorption images of both species following interaction, 
taken for the case of deep confinement (s z = 50) and a velocity v re \ = 2vr. The target's 
velocity distribution is very broad, as expected for a deep ID Mott state, and is centered 
around the velocity of the moving lattice (—2^). The velocity distribution of the probe, 
initially centered around v p = 0, displays a peak of atoms transferred to —4:Vr. A series 
of such scattering spectra (integrated along the transverse y axis) is shown in Fig. [2] (c), 
for a range of relative velocities v re \. For each velocity, the time T move is chosen so that 
the target crystal enters the probe and then moves for 1.4 ms (kept much less than the 
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trap period T z = 2tt/uj z ~ 14 ms to avoid the dispersion of out-coupled atoms to different 
velocities due to propagation in the trapping potential). From the scattered probe spectra 
as in Fig. [2] (b) we observe that, in addition to the line of "transmitted" atoms at v p = 0, 
some probe atoms are out-coupled (i.e. reflected) due to the target crystal. The reflection of 
probe atoms comes from two apparent elastic mechanisms - specular reflection and resonant 
Bragg diffraction. We note that both of these mechanisms are to be distinguished from 
observations of free-space four- wave mixing in two-component mixtures [2Ql EI] , which does 
not persist when one species loses matter- wave coherence. 

The specular contribution is due to the reflection of probe atoms from the potential "step" 
of the target crystal (energy mismatch outside and inside the crystal). This reflection is 
dominant at low velocities and gives way to transmission at larger probe velocities (kinetic 
energies), where the band structure of the target "atomic crystal" is nearly free-particle-like. 
The specular contribution leads to the transfer of probe atoms to v p = —2v re \ and shows 
no resonant structure. In contrast, Bragg reflection from scattering centers of the target 
crystal occurs at values v re i/^R = n and results in the transfer of probe atoms to velocities 
v p = —Vrei — nvR. Such a "locking" to a Bragg resonance is observable through a change 
in the slope of the out-coupled branch or probe atoms. Results of the velocity scan, which 
displays these features, are shown in Fig. [2] (c). 

In Fig. [2] (d) we count the number of probe atoms out-coupled to v p ~ 2v re i ) and plot it 
as a function of v re i/vR. Three resonant peaks in the number of reflected atoms are observed 
for v re i/vR = 1, 2, 3, corresponding to first-, second-, and third-order Bragg processes. In ad- 
dition to these resonances, a significant non-resonant contribution due to specular reflection 
is observed, which as expected decays with increasing velocity v re \. We fit the second-order 
resonance peak with a Gaussian on top of a linear slope [red line in Fig. [2] (d)], and extract 
a 1/y/e- width of 2a v = 0.3vr, reflecting the in-trap velocity width of the probe species. The 
width is a factor of two less than the directly observed TOF velocity width of the probe 
atoms [29J , likely owing to effects of interaction-induced expansion [30] . 

Using our ability to vary the total time T move that the target crystal moves at constant 
velocity, we can also study the temporal buildup of the out-coupled probe population. In 
contrast to optical scattering, this may in the future be used to provide "line-of-sight" 
information in matter- wave scattering experiments. In Fig. [2] (e), we plot the temporal 
buildup of probe population to v p ~ —2v re i, for the cases v re \ = 1.5^ and 2vr. The 

7 



difference between incoherent (specular reflection) and coherent (Bragg) processes should 
lead to differences in growth behavior. Indeed, there are some indications of a more nonlinear 
initial growth for v re \ = 2vr than for v re \ = 1.5^ case. However, since the relatively short 
coherence length of the probe [L ~ h/(ma v ) = 0.8 /im] precludes a fully coherent temporal 
evolution, the observed growth will be mostly determined by the time-dependent overlap of 
the inhomogeneous profiles of the probe and target. 

Finally, we recall the dependence of the static structure factor S(q) on the atomic form 
factor \f{q)\ 2 ) in relation to the three Bragg resonances observed in Fig. [2] (d). In the 
harmonic approximation, the form factor is proportional to exp(— q 2 a 2 /2) ) where a z is the 
harmonic-oscillator length characterising the extent of the atomic distribution on each site. 
While the third-order peak (q = 3) is smaller than the second-order peak (q = 2), the ex- 
pected momentum dependence is partially masked, owing to a fixed probe-target interaction 
time that results in the probe atoms of lower velocity (such as for the first-order peak) in- 
teracting with fewer scattering centers of the target crystal. However, we can directly probe 
the contribution from the atomic form factor, which can be seen to be formally identical to 
the Debye-Waller factor (exp(— q 2 (u 2 )/2)) describing the reduction of elastic scattering due 
to position fluctuations (u 2 ) of scatterers in an ionic crystal pp. The larger a Z) the smaller 
\f{q)\ 2 ) in analogy to the decrease of the Debye- Waller factor with temperature. We point 
out, however, that in our case these fluctuations are not thermal but arise from zero-point 
motion, and may be tuned through o z by varying s z . 

Using the earlier described method of probe acceleration in the longitudinal trapping 
potential (with |1,0) target atoms as before), we can study the dependence of the second- 
order Bragg diffraction amplitude on lattice depth s z , as shown in Fig.[2](f). The plot shows 
an increase of diffraction intensity as the target atoms become more localized, and the data 
is seen to be in good qualitative agreement with the calculated dependence of \f(q)\ 2 on 
s z (using Bloch functions for (fo). Even better agreement can be obtained by taking into 
account ^-dependent heating due to Rayleigh scattering from the laser beams forming the 
longitudinal lattice (see Methods). 

Detecting forced antiferromagnetic ordering 

To further demonstrate that matter- wave diffraction may be used to probe non-trivial struc- 
ture, we show here that it gives access to a clear signature of forced antiferromagnetic 
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ordering [18] in a mixed-spin crystal. As before, we start with a mixture of atoms in one- 
dimensional tubes, with a probe that sees no longitudinal lattice along z here |1,1) 
atoms, with X z <~ 788 nm). However, instead of using only one additional species that 
experiences an attractive lattice along z, our target now consists of two species. The first 
one |r) = |1, — 1) again experiences an attractive lattice potential, whereas for the second 
species \b) = |2, —2) the potential is repulsive. Thus, the |r) atoms will be drawn to the in- 
tensity maxima, while the |6) atoms will be forced to the minima, as illustrated in Fig. [3] (a). 
While we keep the total population of these two species fixed, we vary their relative popu- 
lation, which can be quantified by the net "polarization" of the crystal P = AN/N, with 
AN = N r - N b and N = N r + N b . 

When the target is "spin-polarized" consisting of either all |r) or all |6) atoms, the situa- 
tion is as before, with a lattice constant of d = X z /2 ~ 400 nm. However, when both species 
are present, the new dominant length scale, between atoms of the two species sitting on 
distinct sets of sites, will be given by d! = A z /4 ~ 200 nm. This may be viewed as a crystal 
with the original periodicity and a two- atom basis, or as a new crystal structure having half 
the lattice period, given that the interactions of the probe with the \r) and |6) atoms are ap- 
proximately the same (see Methods). This change results in a different diffraction spectrum 
for the matter wave probes. With the addition of a second species, as scattering centers 
separated by d give way to those of a smaller spacing d', the original first-order diffraction 
peak at v re i/^R = — 1 will become diminished, and entirely disappear if an equal mixture of 
the two species uniformly fills the lattice, with a new first-order peak now occurring at half 
the momentum-space frequency, v re i/^R = — 2. 

To probe the mixed-spin crystal, we move it at a constant velocity of v re i/vR = — 1 with 
respect to the probe atoms, for a fixed interaction time of 1 ms. In Fig. [3] (a), we show 
probe velocity distributions for the three cases P ~ — 1, 0, and 1. As can be seen, the 
number of probe atoms out-coupled to 2vr is much lower for the spin mixture (P = 0) 
than for either of the nearly spin-polarized cases. We note that an appreciable number of 
out-coupled atoms appears even for the balanced spin-mixture, which is most likely due to 
specular reflection as in the previously studied case, while the presence of unoccupied sites 
of either the attractive or repulsive lattices may also cause some Bragg diffraction consistent 
with the original lattice spacing. 

To characterize how the crystal structure changes as the population imbalance is contin- 
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uously tuned, we count the number of probe atoms transferred to a velocity region around 
v p ~ 2vr. As a function of the crystal "polarization", the transferred population shows a 
distinct minimum near P = for a balanced mixture, as shown in Fig. [3] (b). This exam- 
ple readily demonstrates how matter-wave scattering can be used to detect changes to the 
crystal structure of an ultracold lattice gas, and through species-selective scattering may 
eventually be used to detect quantum-magnetic spin-ordered states. Finally, we point out 
that the full tunability of the wavelength for atomic matter waves should provide a dis- 
tinct advantage over optical scattering in certain circumstances, namely in studying density 
structures not formed in optical lattice potentials, such as through self-organization or with 
magnetic trapping potentials. By using probes of small de Broglie wavelength, this could 
also include the study of novel quantum states realized at more easily attainable tempera- 
tures in systems with smaller characteristic length scales. For example, at fixed lattice depth 
s and tunneling-to-interaction ratio t/U (e.g., by control of the scattering length a 5 ), the 
Neel temperature T N) marking the onset of antiferromagnetic ordering, scales with lattice 
spacing d as T N oc d -1 / 2 [3]], and has an even more favorable scaling for fixed a s and variable 
s. 

Conclusion. 

In this work, we have demonstrated that matter-wave diffraction can be used to characterize 
the crystalline structure of strongly correlated atoms in an optical lattice. In the future, 
these techniques may be extended to the characterization of various novel states of ultra- 
cold matter, such as charge- and spin-density waves, magnetically-ordered ground states of 
quantum gas mixtures, and even self- assembled structures such as Tonks-Girardeau gases 
of fermionized bosons, Abrikosov vortex lattices [321 [33], and dipolar crystals [3^] . 

Methods 

Mixture preparation and detection. Starting with an optically-trapped Bose-Einstein con- 
densate of 87 Rb atoms in the |F, mp) = |1, —1) hyperfine ground state, the system is split into an array of 
one-dimensional tubes in the horizontal xy plane. Next, at a magnetic bias field of 1.7 G along z, we create 
hyperfine state mixtures via combinations of microwave Rabi pulses and Landau-Zener sweeps. Different 
combinations of hyperfine states are used to perform particular experiments on interspecies scattering, with 
slight variations in the mixture characteristics. For detection, all species are absorptively imaged on the 
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F = 2 — >> .F' = 3 cycling transition, concurrent with optical pumping from F = 1 — >> F' = 2. 

The initial two-species mixture of probe |p) = |2, —2) and target |t) = |1, 0) atoms, used to study the s z - 
dependent scattering [Fig.[l]; Fig. [2] (f )] , contains a total of (1.6 ±0.2) x 10 5 atoms, with 60% of atoms in the 
target state. For this mixture, a final magnetic bias field of 7.4 G is employed to suppress |1,0) |1,±1) 
spin-changing collisions. A second binary mixture of probe \p) = |2, — 2) and target \t) = |1,— 1) atoms 
is used in conjunction with the moving optical lattice [Fig. [2] (a-e)]. This mixture contains a total of 
(2.8 ± 0.5) x 10 5 atoms, with 33% target atoms. A three-species mixture of one probe species (\p) = |1, 1)) 
and two target species (|2, —2) and |1, —1)) is used to study forced antiferromagnetic ordering [Fig. [3]. This 
mixture contains a total of (1.5 ± 0.2) x 10 5 atoms, with 50% of the atoms in the target. The target, 
consisting of two different species, has a fully tunable spin composition [c.f. Fig. |3](b)]. The intra- and 
interspecies scattering lengths for all the states used are approximately equal to 5.3 nm, i.e. the background 
scattering length for 8T Rb atoms. 

One-dimensional tubes. The array of one-dimensional tubes is created by ramping up two 
orthogonal lattices (along x and along y) within 200 ms to depths of s± =40, via partial retroreflection of 
the ODT's laser beams [35], resulting in a final trapping frequency uj z /2tt ~ 70 Hz along z as determined 
by dipole oscillations. In the harmonic approximation, the transverse oscillation frequency in each tube is 
u±/2tt = 26 kHz, resulting in a spacing to the first allowed transverse excited mode of 2huo± « 14 Er, which 
is greater than the probe kinetic energy T = {v^/vr) 2 Er and all other relevant energy scales (thermal as 
well as all intra- and interspecies interaction energies). For all the cases studied the probe species, while 
ID, is not deeply within the Tonks-Girardeau regime [36], and is characterized by a Lieb-Liniger parameter 
value 7 < 1 [371 EH] 

State-selective lattice potential. The state-selective lattices are formed by interfering two laser 
beams (1/e 2 radius ~ 230 /urn , the same polarization) with tunable wavelength between the 87 Rb Di and 
D2 lines to effect a light-shift cancellation for the probe atoms. In the case of a stationary target of |1,0) 
atoms, the lattice is loaded to a variable depth s z with an s-shaped curve in 75 ms and held for an additional 
5 ms prior to acceleration of the probe atoms. This lattice is made from fully retroreflected laser light of 
wavelength A = 785 nm of a~ polarization. For the case of a moving target of |1, —1) atoms (wavelength and 
polarization as in the stationary case), the lattice is first smoothly loaded in 45 ms to a depth of s z = 10, 
exceeding the critical depth of the ID Mott insulator transition, and then loaded in 5 ms to a depth of 
s z = 50 to freeze the atoms to the sites of the lattice. The lattice is then moved, by introducing a relative 
frequency detuning of Suj between the forward and retroreflected laser beams comprising the z-lattice, via 
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two acousto-optic modulators driven by phase-locked function generators. 

For the case of a spin-mixed crystal with forced antiferromagnetic ordering, the two species are quickly 
loaded into the state-dependent lattice in 2 ms following creation of the target mixture to avoid spatial 
separation in a small magnetic field gradient used to separate the probe. The lattice has modulation depths 
s z = 5 and s z = 12, for the \b) = |2, — 2) (repulsive lattice) and \r) = |1,— 1) (attractive lattice) atoms, 
respectively. This lattice is formed with light of wavelength A ~ 788 nm, of slightly elliptical polarization. 
After loading, this lattice is moved at a fixed velocity as above, with a restriction to relatively low velocities 
to ensure that both target species faithfully follow the optical potential at these modest depths, as observed 
in the velocity distributions of the target species. 

All the lattice depths are calibrated using Kapitza-Dirac atom diffraction [39] , with a systematic uncer- 
tainty of about 3%. For probe atoms, the lattice potential is sufficiently "zeroed" even for greatest available 
optical potentials such that no diffraction is observed for a pulsed-on lattice or for application of a linear po- 
tential gradient while the lattice is present (i.e. no Bloch oscillations). To further "zero" the optical lattice, 
the probe atoms are loaded into ID tubes, and we either minimize excitations of probe atoms undergoing 
dipole oscillations, or in the case of a moving lattice we move with velocity of vl = vr, with no observable 
transfer of probe atoms to non-zero velocities. 

Second-order Bragg peak population vs. s z . The percentage of second-order Bragg-reflected 
probe atoms Ndiff plotted in Fig. [2] (f) is determined from a fit of the TOF spectrum with Gaussian peaks 
for (1) the transmitted probe around v p « v® ; (2) a broad incoherent background centered at v p = ; (3) 
Bragg-reflected atoms at v p = —2vr. N^iff is determined as the amount of Bragg-reflected atoms normalized 
with respect to the total probe population. 

Heating effects. The observed percentage of Bragg diffracted atoms is in qualitative agreement 
with that expected from the form factor \f(q)\ 2 , up to an overall scaling factor [dashed line in Fig. |5](f)]. 
However, the observed signal saturates at large s z while the expected curve does not. We attribute this 
mostly to a heating of the probe species, leading to reduced coherence and an increased spectral width, as 
observed through a linear increase in the probe's TOF velocity-width with increased hold time of a deep 
lattice (s z — 25.5) prior to probe acceleration to v® = 2vr. The Bragg diffraction signal exhibits a roughly 
exponential decay with hold time(l/e-time theat ~ 150 — 200 ms). Assuming that contributions to this 
from heating due to Rayleigh scattering should scale as the time integral over the lattice depth, we expect a 
correction factor of exp(— 0.01s z ) to modify the form factor. A fit to the data with \f(q)\ 2 x ex.p([3s z ) yields 
(3 = —0.015 [value used for the dashed curve in Fig. |5](f)], consistent with the role of such heating. 
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FIG. 1. Interspecies collisions of one-dimensional bosons, (a) Probe atoms of species \p) are 
displaced by a distance Az and then accelerated within a trap potential of trap frequency uo z to a 
final velocity v® = oj z Az. The \p) atoms then interact with the stationary target atoms \t) at the 
trap center, (b) Time-of-flight (TOF) absorption images of the target (top) and probe (bottom) 
species, with velocity components at v p j — and 2vr, due to momentum-exchanging "Newton's 
cradle" (NC) collisions, (c) Illustration of momentum and energy exchange in the band structure 
of the optical lattice, (d) Percentage of target atoms that participate in NC-type collisions. The 
solid curve is an exponential fit with a decay constant of 0.35 Er. The percentage is determined 
from the peak asymmetry in the TOF velocity distribution (after summation over the y direction) 
of the target species, as shown for v® — 2vr in the inset, with atom numbers N± 2 (peaks) and 
TV^ot (total). 
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FIG. 2. Probe scattering from a crystalline target, (a) Target atoms strongly confined to a 
moving state-selective lattice, impinge upon a stationary cloud of probe atoms \p) at a velocity 
vl = 5(jj(\/4:7v), where 8uj is the frequency detuning between the two lattice beams, (b) TOF 
images show the \t) and \p) species after interaction, (c) Probe TOF spectra after interaction 
with a target at s z — 50, at varied relative velocity v re \. Each horizontal line is obtained from a 
TOF image as in (b) by integrating along y and normalizing to the total \p) atom number. The 
solid white lines near v re i/vR ~ 2 illustrate the slope of the clearly visible second-order Bragg 
resonance. The lines are extrapolations of a linear fit to the positions of maximum out-coupling for 
1-8 < Vrel/ V R ^5 2.2. (d) Dependence of reflection from the crystal on v re \. We count the number 
of atoms scattered to v p ~ —2v re i, i.e. within dashed lines shown in (c). Bragg resonances of the 
\t) crystal are observed for v re i/vR = 1,2,3, on top of a background contribution due to specular 
reflection, which decays with increasing v re \. (e) Reflected probe population as a function of the 
time during which the crystal is moved (T move ), for fixed velocities v re l/vR = 1.5 (black data) and 
Vrel/ V R — 2 (red data). The data fall off at long times due to evolution in the trapping potential, 
(f) Percentage of Bragg-diffracted probe atoms (with v® = 2vr) as a function of target confinement, 
s z . The dashed curve is proportional to the form factor |/(2 x 2n/d)\ 2 calculated as a function 
of s z ; the solid curve accounts for additional heating effects (see Methods). The black square is a 
reference measurement without \t) atoms. Error bars represent one-standard deviation statistical 
uncertainties, with thin-line error bars in (e) representing extrapolated errors for individual data 
points. 
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FIG. 3. Detecting forced antiferromagnetic order via matter-wave scattering, (a) Velocity dis- 
tributions of scattered probe atoms following interaction with a moving (at v re \ — —vr) crystal 
of atoms in a tunable spin-mixture of red \r) and blue \b) states, attracted to and repelled from 
intensity maxima of a state-dependent lattice, respectively. Shown are the cases of a crystal of 
nearly all \b) atoms (P ~ —1), nearly all \r) atoms (P ~ 1), and a balanced mixture of the two 
states (P ^ 0). (b) The number of Bragg-diffracted probe atoms [shaded regions in (a)], as a 
function of the spin population imbalance of the target crystal, P = AN/N. A distinct minimum 
is observed for P « 0, relating to the case of a balanced spin- mixture crystal with maximal forced 
antiferromagnetic order at spacing d! — A z /4. 
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